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ABSTRACT 

Context. The '^C/'^'C isotope ratio in the interstellar medium (ISM), and its evolution with time, is an important tracer of stellar 
yields. Spatial variations of this ratio can be used to study mixing in the ISM. 

Aims. We want to determine this ratio and its spatial variations in the local ISM from CH+ absorption lines in the optical towards 

early-type stars. The aim is to determine the average value for the local ISM and study possible spatial variations. 

Methods. We observed a large number of early-type stars with Feros to extend the sample of suitable target stars for CH+ isotope 

studies. The best suited targets were observed with Uves with higher signal-to-noise ratio and spectral resolution to determine the 

isotope ratio from the interstellar CH^ lines. This study significantly expands the number of '""CH^ detections. 

Results. We find an average ratio of (R) = 76.27 ± 1.94 or, for / = 1 /R, </> = (120.46 ± 3.02) x 10"-*. The scatter in / is 6.3 x cr«/)). 

This findings strengthens the case for chemical inhomogeneity in the local ISM, with important implications for the mixing in the 

ISM. Given the large scatter, the present-day value in the ISM is not significantly larger than the solar value, which corresponds to the 

local value 4.5 Gyr ago. 
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1. Introduction 

, The '^C/'^C ratio is a cornerstone of models of the nuclear his- 
' tory of our galactic interstellar matter (ISM), that is, the Galactic 

■ Chemical Evolution. Low and intermediate mass stars increase 
. (during the AGB phase) the '^C/'^C ratio to ~ 300, while stars 

■ more massive stars (undergoing CNO processing) increase the 
' ratio to '^C/'^^C~3. A measurement of the '^C/'^^C ratio in the 

ISM thus provides important information about the relative im- 
portance of stars of different mass to the chemical enrichment of 
the interstellar matter. 

The '^C/'-'C ratio in the solar system is 89 (Anders & 
Grevesse 1989). Models of chemical evolution predict a decrease 
in the '^C/'^^C ratio with time for a given galacto-centric dis- 
tance, and a decrease with galacto-centric distance in the Galaxy 
(Prantzos et al. 1996). Such models assume that mixing in the 
ISM is complete and restricted to material at a given galacto- 
centric distance (azimuthal mixing). A comparison of results be- 
tween our Galaxy and other galaxies will give important data 
for the nuclear processing history (see, e.g., Tosi 2000; Prantzos 
2001) 

The '^C/'^C ratio has been determined mostly from radio- 
astronomical measurements of CO and H2CO (see e.g., Wilson 
& Rood 1994) and millimeter observations of CN (Savage et al. 
2002; Milam et al. 2005). 
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These results may be affected by interstellar chemistry in 
different ways, however: Chemical fractionation (which results 
from different zero-point energies of the isotopomers) enriches 
molecules in '^C (Watson et al. 1976) and thus lowers the ra- 
tio, selective dissociation destroys the rarer species more and 
thus raises the ratio. The latter effect is due to self-shielding of 
the more abundant isotopemer, which therefore survives photo- 
destruction better than the less abundant counterpart. 

The molecule CH^ is thought to be produced only in hot- 
ter parts of photo-dissociation regions and is a rare species, so 
is not affected by either chemical fractionation or selective dis- 
sociation. In addition, the optical CH^ lines have typically only 
small saturation effects. This molecule, which can be observed 
in optical absorption lines from the ground, is therefore consid- 
ered to provide the most secure '^C/'-'C ratios (see e.g. Watson 
et al. 1976). We caution, however, that the formation process of 
CH^ is probably not completely understood (Gredel et al. 1993). 

Two systems of CH^ lines are strong enough to be easily ob- 
served from the ground. In addition, the absorption lines from 
'^CH+ at 4232 and 3957 A are sufficiendy shifted from their 
'^CH^ component to allow unblended measurements, provided 
the velocity profiles in the intervening cloud are narrow enough. 
If both are measured simultaneously, we can make use of the 
fact that the '^CH+ line is blue-shifted by 0.26 A relative to the 
12CH+ at 4232 A but red-shifted by 0.44 A relative to the ^^CH^ 
line at 3957 A. This way, a chance superposition with another 
cloud along the line-of-sight can be excluded. In addition, the 
two simultaneous measurements help to eliminate systematic ef- 
fects. 
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This property has been used e.g. by Stahl & Wilson (1992) to 
check that there is no accidental overlap of a weak velocity com- 
ponent of i^CH"" at -18.8 km sec'^ from the deepest ^^CH"" ab- 
sorption line with the ^^CH"^ line. Nearly all previous data were 
taken at 4232 A only. 

These measurements in the optical require, however, obser- 
vations with high spectral resolution and signal/noise ratio. Most 
of the published CH"^-measurements have therefore been ob- 
tained with very bright background stars at relatively small dis- 
tances from the Sun. 

The availability of efficient high-resolutions spectrograph at 
8m-class telescopes in recent years has improved the situation 
considerably (Casassus et al. 2005, in the following paper 1). 
The observations are hampered, however, by the scarcity of suit- 
able pubUshed candidate stars with sufficiently strong and nar- 
row CH"^ absorption lines. In this paper we describe new obser- 
vations which have been obtained at targets specifically selected 
for this purpose. 

The outline of the paper is as follows: Section 2 describes 
our observations. Section 3 gives details on our procedure to 
measure the isotopic ratio R. The results are discussed in sec- 
tion 4, both the average properties of the sample and notes about 
individual targets. Section 5 discusses the results and concludes 
the paper. 

2. Observations and reduction 

2.1. Preselection with Feros 

Many of the pubUshed CH"^-measurements are hampered either 
by too broad or complex line profiles of the interstellar Unes. 
The isotope shift is very small, so that for most sight-lines (with 
relatively broad absorption lines) the isotopic lines partly over- 
lap, which introduces considerable errors to the isotope ratio. 
Interstellar lines with several components are even more prob- 
lematic. The best lines of sight have sharp, single absorption 
Unes. 

The stars with published CH^ -profiles that are suitable for 
the determination of the '^C/'^C isotope ratio are exhausted by 
previous studies. We therefore undertook a dedicated search for 
additional targets, which might be suitable for studies of the iso- 
tope ratio. 

For this purpose we used the Feros echelle spectrograph 
at the ESO/MPl 2.2m telescope at ESO, La Silla. Feros has 
a sufficiently high spectral resolution (48 000) to examine the 
CH"^ line profiles. In three nights in January 2006 we observed 
about 80 early-type stars with a B magnitude brighter than 8.0 
selected from various lists (Humphreys 1978; Maiz-Apellaniz 
etal. 2004). 

The exposure times for the Feros spectra were relatively 
short (typically a few minutes), just long enough to give a rea- 
sonable S/N-ratio near the ^4232, 3957 Unes of CH+. The ^^CH+ 
component is not detected in these spectra. The Feros data were 
pipeline-reduced directly at the telescope. From these spectra we 
selected the targets which were best suited for follow-up obser- 
vations with UvES. The main selection criteria were: 

- sharp and symmetric CH^ -lines 

- high equivalent width and 

- brightness of the target in the B band 

The sharpest and deepest absorptions are needed to separate 
"CH-'A4232 from i2CH-'A4232. We first selected by inspection 
the narrowest and single-component absorptions, with line pro- 
files close to a single Gaussian. We then ordered individual lines 



of sight according to two quantitative indicators of sharpness. 
The first indicator is the signal-to-noise ratio on the absorption 
Une, r = D/o- = VF^/A = W^lO' ^-^^/A, where D is the 
depth of the line, Wa = ( ^/7T/(2 V2 In 2))DA/Fc its equivalent 
width, and A its full width at half maximum, cr oc -^7^ is the 
noise, and B is the B-band magnitude of the background star. 
The second indicator, introduced to weight the final ranking to- 
wards the sharpest lines, was simply = y/A. The final list of 
stars selected and later observed with Uves is given in Table 1 . 

Since the Feros spectra cover a large spectral range (3600 
- 9200 A) we also obtained information about many other in- 
terstellar lines. A paper describing the fuU Feros data set is in 
preparation. 

2.2. High S/N observations with Uves 

Observations of the selected stars with high S/N were obtained 
with the Uves echelle spectrograph at the VLT unit telescope 
Kueyen at Cerro Paranal, Chile, in three nights between January 
8/9 and January 10/11, 2006. Uves allows measuring both line 
systems in the same detector setting, thus providing a means to 
correct for line blending if it is apparent in the spectra. 

The observations are difficult because of the need for both 
very high S/N and very high spectral resolution to acquire the 
profile of faint '^CH'*", against a very bright stellar continuum 
and close to the much stronger '^CH"^ component. We need S/N 
ratios of at least ~10 on '^CH"^ for accurate profile fitting. 

The project requires the highest possible spectral resolution. 
We therefore used an image slicer to minimize flux losses. In ad- 
dition, the image slicer distributes the light along the slit, which 
improves flat-fielding and aUows longer integration times before 
saturation occurs. Uves slicer #2 was used, which reformats an 
entrance opening of I'.'E x 2'.'2 to a slit of 0'.'44 width and a 
length of 7''9, which is imaged on the spectrograph entrance sUt 
of 0'.'45 width. The spectral resolution in this configuration is 
A/AA = 75 000. The central wavelength was set to 4 370 A, which 
gives a spectral coverage from 3 730 to 5 000 A. During our ob- 
servations the seeing was sometimes too good to fiU the entrance 
aperture of the images sUcer. Therefore, for some spectra, only 
one or two sUces contained most of the signal, which unfortu- 
nately decreases the gain of the image slicer. 

Accurate flat-fielding is also important. Therefore a large 
number of flat-fields (150) was obtained during day-time dis- 
tributed along the observing run. 

In addition, two rapidly rotating unreddened early-type stars 
(HD 10144 = ffEri, spectral type B3Ve and HD 108248 = a'Cru, 
spectral type B0.5IV) were observed in order to check for pos- 
sible faint telluric features. These stars have been observed with 
exactly the same instrumental settings as the target stars and with 
similar S/N-ratio. We confirm there are no detectable telluric fea- 
tures under the CH"^ absorption. 

The brighter targets where exposed until about 50% of the 
maximum level allowed by the CCD detector was reached. This 
to typically a few minutes for our targets. Series of up to 50 
exposures per target were obtained to build up the required S/N- 
ratio. Targets which turned out to show too complex line profiles 
were dropped from the target list. The observations are summa- 
rized in Table 1, which also lists the mean airmass of the target 
objects in the different nights. At the higher airmass values, the 
range in airmass differed from the mean by about ±0.1 during 
the observations, and less at smaller airmass. 

We used the UvES-pipeline software for the reduction of the 
spectra. In order to maximize the S/N-ratio of the extracted 
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spectra, all flat-fields obtained during the run were averaged. 
A mean ThAr- spectrum obtained during day time was used for 
the calibration of all spectra obtained in one night. After back- 
ground subtraction and flat-fielding, the spectra were extracted 
as 2D-spectra using a long extraction slit covering all slices. All 
spectral orders of these spectra were then merged in one 2D- 
spectrum, keeping the information along the slit. Finally, the 
spectra were collapsed to a ID-spectrum by averaging along the 
slit and thereby averaging the flux from all slices. The region 
around the CH^/13957 A-line was influenced by a detector blem- 
ish. This defect only covered two CCD lines along the slit. In or- 
der to remove this defect, we also extracted spectra excluding the 
two afifected CCD lines. The observations of of each target were 
then averaged to a mean spectrum. Wavelengths are reported in 
air and referred to the solar barycenter. 



Table 1. Summary of the observations. The table columns contain the 
HD number, spectral type, B magnitude, average airmass, total inte- 
grated exposure time and tlie number of exposures for each target. 



Target 


Spec. 


B Airmass 


Exp. 
LsccJ 


Number 


HD10144 


B3Ve 


0.3 


1.36 


63 


90 


HD23480 


B6IVe 


4.11 


1.52 


1100 


55 


HD35149 


BIV 


4.85 


1.16 


3920 


98 


HD37903 


B1.5V 


7.91 


1.13 


5760 


24 


HD52266 


09V 


7.22 


1.15 


606 


3 


HD52382 


Bllb 


6.64 


1.05 


4800 


40 


HD53974 


B0.5IV 


5.44 


1.03 


3020 


60 


HD58510 


Bllab 


6.91 


1.04 


5000 


25 


HD73882 


09III 


7.58 


1.05 


6257 


27 


HD75149 


B3Ia 


5.68 


1.09 


5400 


90 


HD76341 


Bl/B2Ib 


7.39 


1.06 


6000 


25 


HD91452 


BOIII 


7.66 


1.30 


5760 


24 


HD92964 


B2.5Iae 


5.59 


1.22 


3900 


65 


HD108248 


B0.5IV 


1.32 


1.31 


20 


20 



3. Model line profiles 

We model the CH^ opacity in velocity space, as in paper 1. But 
in this work we force the same opacity profile t(v) on both over- 
tones. In paper I we preferred to fit the two overtones separately, 
and thus obtain independent measurements with which to assess 
the role of systematics. Here we incorporate the prior knowledge 
that both CH^/l/l3957,4232 overtones share a common opacity 
profiles, which is true because they both stem from the same de- 
excited energy level (the ground state, see Casassus et al. 2005, 
Sec.l) 

Requiring a common velocity profile highUghts the limits 
of our knowledge of the underlying continua on the wavelength 
scales of the CH"^ Une widths. The continua can be affected by 
detector glitches, or atmospheric features. In the absence of in- 
strumental artifacts or stellar features, the two profiles should 
be the same. Any difference in the fit residuals should highlight 
problems with the continuum definition. 

In this paper we investigate how a simultaneous fit to both 
overtones helps constrain the CH^ opacity profiles, decomposed 
in a set of Gaussians in velocity. This approach was attempted 
in paper I to resolve cases when the '^CH"^/14232 transition was 
blended with the main isotope. But in paper I the quahty of the 
simultaneous fits was not satisfactory. We now explain this prob- 
lem in paper 1 as being due to our overlooking the possibihty of 



small shifts in the observed relative wavelengths of the two over- 
tones, AAo. Allowing for AAo produces satisfactory simultaneous 
fits to both overtones. Such wavelengths shifts are of the order of 
the uncertainties in the absolute wavelength calibration of each 
echelle of a few mA. From the fits to the 1 1 fines of sight ac- 
quired in 2006, we find that the maximum range of values for 
A/lo is 3.2 10"^ A, and its Icr scatter is 9.5 10"''A, while the 
resolution element is 5.3 10"^ A. 

The isotope ratio information in cases where the rare isotope 
is blended at /14232 cannot be recovered without the /13957 in- 
formation. The uncertainties in the baseline definition make it 
impossible, using a single transition, to distinguish the rare iso- 
tope from a detector gUtch. But AA22>2 blends are still useful 
as a consistency checks. We incorporate the uncertainties in the 
baseline definition under /14232 through the inclusion of extra 
free-parameters. The possibility of detector glitches or stellar or 
atmospheric features is modeled with a correction to the baseline 
under the interstellar absorption. We explored the possibilities 
of two broad types of low-level baseline corrections: 1- additive 
corrections, as in low-level atmospheric, stellar or CH^ features, 
and 2- multiplicative corrections, as in flat fielding features. The 
low-level baselines were modeled using either an expansion in 
Legendre polynomials or as segmented fines. But the Legendre 
expansion required high orders, and thus many free-parameters, 
because of the condition of fixed endpoints imposed by continu- 
ity with the baseline neighboring the interstellar absorption. We 
obtained best results with an additive baseline correction, com- 
posed of a segmented line defined by the ampUtude of spikes 
distributed uniformly in wavelengths. 

We caution that the present technique used to measure the 
CH"^ isotopic ratio suffers from the same uncertainties inher- 
ent to any observation of absorption lines. The inferred equiv- 
alent widths depend on the exact placement of the continuum 
level. The uncertain definition of the baseline is a source of 
systematic uncertainty. However we take a conservative value 
for the thermal noise, with the hope of accounting for at least 
small variations in the baseline definition. Our noise estimate 
corresponds to the standard deviation of the residuals in the 
whole 1.5 A windows around CH'^/l/13957,4232, after baseline 
and profile fitting, including detector gUtches and any telluric 
features. Provided the baselines are realistic enough, our proce- 
dure should indeed allow for small variations among plausible 
baseUnes. 

The model we use to define the underlying continua and to 
parameterize t(v) is otherwise the same as in paper 1. For clarity 
we now itemize our new procedure: 

1. extract a 1.5 A spectrum F{X) centered on the CH"^ line, and 
fit a baseline with a 4 - 7* order Legendre polynomial; in- 
clude weights to improve the quality of the baselines near the 
features of interest; store the RMS dispersion of the residuals 
as the spectrum's noise, crp; 

2. define a model spectrum with 

FmW) = ^'c,W)exp(-T(a)), (1) 

where the line absorption opacity t is a superposition of Wg = 
1 to 4 Gaussians common to each isotope, with ^^CH+ com- 
ponents sharing the parameters of the '^CH"^ components, 
except their centroids are translated by a velocity shift Aviso, 
and their amplitude are scaled by a factor / =^^CH"^/^-^CH'^, 
kept as a free parameter: 
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T(A) = YjT°exp{-0.5(A-A°f/cTi) + 

/=1 

no 

^/T°exp(-0.5(i-^°(l+Av„Jc))2/^°f); (2) 



3. The wavelength shift between each isotope is taken as a free 
parameter, allowed to vary over a resolution element. 

4. In cases of blends or suspected glitches we include an ad- 
ditive correction to the baseline (details are given in the de- 
scription of results from individual lines of sights). The cor- 
rections applied to the continuum level derived from Step 1, 
are parameterized using a segmented line defined by «spikes 
points distributed in equal wavelength intervals under the 
whole CH+ profile. 

5. Since blends always occur under /14232, we artificially in- 
crease the noise under the fundamental tone by a factor of 
10, thereby allowing the otherwise noisier overtone to dom- 
inate the optimization. 

6. Optimize the model parameters by minimizing - 
2^(F(/lj)-Fni(/lj))^/crp, using the noise from Step 1. The op- 
timization was carried with the PDL::Minuit package (based 
on the Minuit optimization library from CERN) 

7. Update the noise of the stellar spectrum by replacing crp with 
the root mean square (RMS) dispersion of the residuals (the 
difference between the observed and model spectra over the 
whole 1.5 A window). 

8. Estimate the uncertainty in model parameters by two meth- 
ods: 

(a) search parameter space for the Ax^ = 1 contour defining 
the region enclosing 68.3% confidence level. 

(b) estimate 1 cr uncertainties from the curvature matrix of 
(i.e. approximate to normal errors). 



4. Results of the fits 

4.1. Average properties of the dataset 

In the following we discuss the results of our fits for all sight- 
lines. In order to increase the number of sight-lines, we include, 
in addition to the targets observed in Jan. 2006 also the best 
targets from paper I: f Oph, HD 110432 and HD 152235. These 
stars have been observed with the same instrument setting but 
were analyzed differently. In order not to introduce any bias in 
our results, we applied the fitting method described above also 
for these stars. The results for R and / differ slightly from the 
results in paper I, but are consistent within the errors. 

The results for the fits of all targets are given in Table 2 and 
are plotted in Fig. 1. 

In paper I we summarized our measurements in terms of 
R =i2q13(- jg preferable to use / = l/R since for / the 

noisiest factor occurs in the numerator. We combine all / mea- 
surements in a weighted average (/) by weighting each line of 
sight by 1/cr^, where cr is the root-mean-square uncertainty on 
/. We find if) = (120.46 + 3.02) x lO"'* . The weighted scatter 
of / values is 6.3 x cr((/)), which means that the scatter is truly a 
property of the distribution, and not merely a result of statistical 
fluctuations about a single mean value. Indeed, applying the 
test for the null hypothesis that all measurements derive from a 
single value, we get that the reduced is 3.02 for 14 data points 
and 1 free parameter (i.e. (/)). Such ax^ value allows us to state 
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Fig. 1. A plot of individual / values for the lines of sight studied in this 
work. The galactic longitude of the targets is also indicated. The thick 
lines are ±lcr uncertainties (2 cr in total length). The hatched box is 
the weighted average of the measurements in the figure, and should be 
representative of the local ISM. The height of the hatched rectangle is 
±1 cr. This value is </> = (120.46 ± 3.02) x lO-^ The corresponding 
average for R values is {K) = 76.27 ± 1.94. 



that the observed distribution of / values is a real scatter with 
essentially 100% confidence. 

The corresponding average for/? values is (/?) - 76.27+1.94, 
and can be compared to the value qIR- 1%.21 ± 1.84 published 
in Casassus et al. (2005), calculated from a smaller sample. 

Contrary to /, the isotope ratio R - l/f is not normally 
distributed, so that the uncertainties on R are not directly related 
to confidence levels. Therefore the ISM statistics over the sample 
of sight-lines must be taken on /. Unfortunately in paper I we 
extracted statistics on R values, rather than / values. Here we 
have corrected this mistake. The conclusions from paper I still 
hold, however The positive skewness in R becomes negligible 
for smaller uncertainties and larger R. In the paper I data the ISM 
scatter in R turns out to be 7% smaller than that in / relative to 
the uncertainties on their weighted average. The variations in / 
towards HD 1 10432 and HD 152235 are 3.4cr (rather than 4.3o- 
in R), and are 3.5cr towards HD 1 10432 and HD 170740 (rather 
than 3.2cr). 

4.2. Notes on individual objects 

4.2.1. HD 23480 

HD 23480 (= 23 Tau = Merope) is a well-known member of the 
Pleiades cluster and the only target in our list which has a pub- 
lished i^c/i^c-ratio (41+9, Hawkins & Jura 1987). We observed 
it, despite its northern location, as a reference target. Because of 
its position the sky, the airmass during our observations was rel- 
atively high. Our result is in disagreement with Hawkins & Jura 
(41 + 9, 1987): We do not confirm their low value of R toward 
HD 23480. The Uves CH^ spectrum (see Fig. 2) appears to be 
free of blends or artifacts. 

4.2.2. HD35149 

A detailed study of the interstellar medium towards HD 35149 
(- 23 Ori) based on high resolution UV spectra has been pub- 
Ushed by Welty et al. (1999). Adamkovics et al. (2003) searched 
the star for rotationally resolved C3, but did not detect this 
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Table 2. Summary of the results of the fits. / is given in units of 10"''. n„ is the number of Gaussian components and W,|(4232 A) and W,|(3957 
A) are the equivalent widths of the respective '^CH^ lines, obtained by integrating the fit solutions. Numbers in parentheses are uncertainties in 
the last digit. The corresponding values for the '^CH^ components can be obtained by multiplying the values for the '^CH^ component with /. 



Target 


R ± cr(R) 


/ = 


tcr(/) 


"g 


W,,(4232 A) 
[mA] 


W,,(3957 A) 
[mA] 


HD23480 


96.5±16.3 


103 


6+17.5 


1 


14.94 (2) 


8.82 (3) 


HD35149 


49.4± 8.2 


202 


2+33.5 


2 


9.48 (3) 


5.29 (3) 


HD37903 


97.5±27.2 


102 


5+28.6 


2 


10.05 (3) 


5.80 (4) 


HD52382 


59.7±13.5 


167 


6+37.8 


2 


19.20 (3) 


11.01 (4) 


HD53974 


91.5±18.8 


109 


2+22.5 


2 


11.76 (2) 


6.68 (3) 


HD58510 


43.7±10.8 


228 


9+56.5 


3 


10.30 (3) 


5.75 (4) 


HD73882 


79.7±14.5 


125 


5+22.9 


2 


17.45 (3) 


9.86 (4) 


HD75149 


102.7±21.0 


97 


3+19.9 


2 


10.22 (2) 


5.73 (3) 


HD76341 


56.2± 8.1 


177 


8+25.6 


2 


38.43 (3) 


23.10(4) 


HD91452 


76.0+23.3 


131 


5+40.2 


2 


8.75 (3) 


4.86 (4) 


HD92964 


146.7±39.3 


68 


2+18.3 


2 


8.06 (2) 


4.46 (2) 


<roph 


80.9± 3.0 


123 


5+ 4.5 


2 


23.42 (1) 


13.82(1) 


HDl 10432 


69.1+ 4.4 


144 


8+ 9.2 


2 


13.74 (2) 


7.88 (2) 


HD152235 


92.2± 4.9 


108 


5+ 5.7 


3 


41.78 (3) 


24.88 (2) 




I , , , L I ^ ill 1 u M I I I , I ^ Uj_i u , , , ^ 1 1 , , , I 

4232 4232.5 4233 4233.5 3957.5 3958 3958.5 3959 



Fig. 2. Spectra of HD 23480. Wavelengths are in A, and flux densities in arbitrary units. The individual Gaussian components comprising the fit, 
prior to folding with the instrumental response, are shown only on '^CH+ as light gray lines, but are omitted from '^^CH+ for clarity. The vertical 
dashed lines mark the line centroids, at the average of the Gaussian centroids weighted by their equivalent widths. The units of the j-axis are 
arbitrary, and are scaled so that the median of the object spectrum is unity. Also shown is a magnified version of the object spectrum, by the factor 
indicated on the figure, and offset for clarity. The residual spectrum is shown on top of the object spectrum, and is also magnified by the factor 
indicated on the figure. The noise used to compute the significance of the fits is labeled on the top right. The noise is the RMS dispersion of the 
residual spectrum, and is in the same units as the y-axis. The height of the shaded rectangles on top of the x-axis indicates the relative weights used 
in the baseline definition. 

molecule towards HD 35149. Measurements of diffuse interstel- 
lar bands towards the star have been published by Thorburn et al. 
(2003). The Uves CH^ spectrum (see Fig. 3) appears to be free 
of blends or artifacts. 

4.2.3. HD 37903 

HD 37903 is the illuminating star of the reflection nebula 
NGC 2023. A detailed study of the interstellar medium towards 
HD 37903, based on ORFEUS UV data, has been pubUshed by 
Lee et al. (2002). This star has previously been observed in the 
CH+-line by Gredel et al. (1993). The Uves CH+ spectrum (see 
Fig. 4) appears to be free of blends or artifacts. 

4.2.4. HD 52266 

HD 52266 is considered a field star, but is projected close to CMa 
OB 1 . de Wit et al. (2004) found evidence for a cluster near the 
star. No detailed study of the interstellar medium in this direction 
has been published so far. 



This spectrum is very shallow because we interrupted the ob- 
servation when finding out at the telescope that the line of sight 
to HD 52266 is very structured, and useless for the determina- 
tion of an isotope ratio. The multiple velocity components are 
nonetheless faint enough to have escaped detection in our Feros 
preselection. The Uves CH^ spectrum is shown in Fig. 5. 

4.2.5. HD 52382 

This star is a member of the CMa OB 1 association and has pre- 
viously been observed in the CH^-line by Gredel (1997). The 
UvES CH^ spectrum (see Fig. 6) appears to be free of blends or 
artifacts. 

4.2.6. HD 53974 

This star is considered a field Be star by Zorec et al. (2005). 
No detailed study of the interstellar medium in this direction has 
been published so far. The Uves CH^ spectrum (see Fig. 7) ap- 
pears to be free of blends or artifacts. 
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4.2.7. HD 58510 



4.2.8. HD 73882 



HD 58510 is a little studied Blab supergiant. No detailed study 
of the interstellar medium in this direction has been published so 
far. 

Although the fundamental '^CH^ tone is only marginally 
blended, the baseline at /14232 shows a break right under '^CH+, 
which is very difficult to account for. The Uves CH^ spectrum is 
shown in Fig. 8. 



This member of Vela OB 1 has previously been observed in the 
CH^-line by Gredel et al. (1993). The star is located in the Gum 
nebula. HD 73882 has been observed in the far ultraviolet with 
FUSE. These observations have been used to study e.g. the inter- 
stellar deuterium abundance (Ferlet et al. 2000) and the interstel- 
lar extinction towards HD 73882 (Sofia et al. 2005). FUSE spec- 
tra of the star have also been published by Snow et al. (2002) and 
Rachford et al. (2002). lUE and FUSE data have been used by 
Sonnentrucker et al. (2007) to determine a '^CO/'^CO column 
density ratio of 25 towards HD 73882, significantly below our 
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value for ^^C/^^C. Since Sonnentrucker et al. (2007) find partic- 
ularly low '^CO/'^CO-ratios in lines of sight which sample cool 
and dense gas, they interpret this finding as evidence for chemi- 
cal fractionation of CO in cool and dense regions. The Uves CH^ 
spectrum (see Fig. 9) appears to be free of blends or artifacts. 

4.2.9. HD 75149 

This star is a member of the Vela OBI association and has pre- 
viously been observed in the CH^-line by Gredel (1997). Gredel 
et al. (2002) and Gredel (1999) also observed the star in the in- 
terstellar CN and C2 lines, respectively. Observations of inter- 



stellar Na I and Ca 11 have been published by Cha & Sembach 
(2000). The Uves CH^ spectrum (see Fig. 10) appears to be free 
of blends or artifacts. 

4.2.10. HD 76341 

Although this star is not in the list of Humphreys (1978), it is 
considered by Reed (2000) as a probable member of the OB 
association Vela OBI. Nichols & Slavin (2004) searched for 
unusually strong excited interstellar Ci in stars seen through 
the Vela supernova remnant, indicating shock excitation, and 
found several stars with unusually strong Ci. Both HD 35149 
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and HD 76341 are included in their sample, but do not show un- 
usually strong C I. 

Although relatively faint, HD 76341 gives us the best mea- 
surement of '^CH^/'^CH^. It has a very sharp and relatively 
strong CH^-line. 

The fundamental '■'CH^ absorption at /14232 appears to be 
blended with the main isotope, as well as affected by baseline un- 
certainties (atmospheric, and instrumental detector glitches). By 
artificially changing the noise under /13957 and /14232, we can 
test the effect of imposing the R value from each overtone. We 
find it is impossible to fit the '^CH^/14232 simultaneously with 
'2CH+/14232 without changing the basehne under '^CH+/i4232 
by deviations from the Legendre polynomial fit (step 1) compa- 
rable to the '^CH^ line depth. On the other hand, if we artifi- 
cially decrease the noise under the overtone '^CH^ absorption 
at /13957, then it is possible to also fit '^CH^ in both overtones. 
Fig. 1 1 shows the result of imposing the '^CH^/13957 profile. It 
is clear that the observed spectral points in '^CH^/14232 do not 
share the same profile as the other 3 lines, '^CH^/14232 is there- 
fore affected by either atmospheric features or detector glitches, 
of the type seen at A = 3858.6 A. The uncertainties on R are 
those obtained by varying only R and keeping fixed all other pa- 
rameters. The UvES CH^ spectrum is shown in Fig. 11. 



away O-type star. No detailed study of the interstellar medium in 
this direction has been published so far The Uves CH^ spectrum 
(see Fig. 12) appears to be free of blends or artifacts. 

4.2.12. HD 92964 

Member of the OB association Car OBI. No detailed study of 
the interstellar medium in this direction has been published so 
far The Uves CH^ spectrum (see Fig. 13) appears to be free of 
detectable blends or artifacts. 

4.2.13. ^Oph 

Because of its brightness and clean interstellar line profiles, 
^ Oph is a prototype star for studies of the interstellar medium, 
also for the C'^C'^ ratio. The Uves CH^ spectrum is shown in 
Fig. 14. For detailed discussion see paper I. 

4.2.14. HD 110432 

This star is also discussed in paper I. The sight-line to 
HD 110432 is one of the cleanest sight-lines studied in paper I. 
The Uves CH^ spectrum is shown in Fig. 15. 



4.2.11. HD91452 

de Wit et al. (2004) unsuccessfully searched for a cluster near 
HD 91452. According to de Wit et al. (2005), HD 91452 is a run- 



4.2.15. HD 152235 

This star is another star from paper! We treated this case as a 
blend, since the '^CH^ profile just overlaps with '^CH^ in the 
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fundamental tone. We exaggerated the noise for /14232 by a fac- 
tor of 3, instead of 10 as in the other blended cases. There is 
possible telluric absorption red-wards of '^CH+/i3957. The Uves 
CH+ spectrum is shown Fig. 16. 



5. Discussion and conclusions 

Our systematic search for stars with clean line profiles in CH^ 
and large equivalent widths, which are suitable for isotope stud- 
ies, has resulted in a much increased sample of reliable '^C/'^C 
isotope ratios from optical observations. 



The VLT-UvES spectra of CH^ absorption in the lines of 
sight towards 14 stars have allowed us to measure the carbon 
isotopic ratio R -^^C/^^C with unprecedented accuracy. The sig- 
nificant extension of the sample of Casassus et al. (2005) allows 
us to strengthen their conclusion that there is truly a significant 
scatter in the local ISM. Averaging our measurements for the 
14 lines of sight gives a value representative of the local ISM: 
{R} - 76.27 + 1.94. This mean ratio is in excellent agreement 
with the local ratio determined from radio-astronomical mea- 
surements of CO and H2CO (76±7, Wilson & Rood 1994), and 
also in good agreement with more recent determinations from 
CN (si68, Milam et al. 2005). This agreement of results obtained 
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Fig. 16. Same as Fig. 2, but for HD 152235. 



from different molecules suggests that chemical fractionation 
and selective dissociation do not play a major role for the CO, 
H2CO and CN molecules. The results of Sonnentrucker et al. 
(2007) suggest, however, that CO may be affected by chemical 
fractionation in cool and dense regions. 

The weighted scatter of / values is 6.3 x cr((/)) which di- 
rectly shows that the '^C/'-'C isotope ratio is variable on the rel- 
atively small spatial scales which we sample. Considering this 
intrinsic scatter of the '^C/'^C ratio in the ISM, the solar car- 



bon isotope ratio of 89 is undistinguishable from the present-day 
value, even after 4.5 Gyr of galactic evolution. This statement is 
not affected by instrumental uncertainties, although a sampling 
bias could distort our ISM value. Of course, this conclusion is 
based on the assumption that the solar system ratio is represen- 
tative for the galacto-centric distance Rq = 8.5 kpc, 4.5 Gyr ago. 
Since we have shown that the interstellar '^C/'^C isotope ra- 
tio is spatially variable now, it is not clear if the solar value is 
characteristic for the local '^C/'^C ratio, 4.5 Gyr ago. However, 
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there is also no strong evidence that the solar value is untypical: 
Although it has been claimed that the sun is overabundant, espe- 
cially in carbon and oxygen, compared to the local ISM, the solar 
carbon and oxygen abundances have been revised downwards. 
Recent papers (e.g. AUende Prieto 2007) find good agreement 
for the carbon and oxygen abundance in the sun, the local ISM 
and young stars. Comets show a '^C/^^C ratio consistent with the 
photospheric solar value (e.g. Manfroid et al. 2005). Therefore 
it seems also unlikely that the '^C/'-'C ratio was modified from 
the value in the solar nebula by fractionation of carbon isotopes 
during the formation of the sun. 

Therefore we conclude that there is no evidence that the solar 
'^C/'^C ratio is untypical for the solar neighborhood 4.5 Gyr 
ago. It is plausible however, that the average value for the solar 
neighborhood was somewhat different from the solar value by a 
similar amount as the present-day variations. A small decrease 
of the average ^^C/'^C ratio in the solar neighborhood during the 
lifetime of the sun is therefore certainly possible. 

The ^^C/^^C ratio increases with galacto-centric distance 
(Wilson & Rood 1994) by about 7.5±1.9 per kpc. Therefore, a 
comparison of the present-day local '^C/'^C ratio with the solar 
photospheric value should also take into account that the sun was 
born at a different galacto-centric distance (Wielen & Wilson 
1997). Since the sun has, according to Wielen et al. (1996), mi- 
grated from its birth-place in the irmer part of the Galaxy out- 
wards by 1.9 ± 0.9 kpc during its lifetime, this effect decreases 
the expected change in the ^^C/^^C ratio significantly. 

Early models for the ^^C/^^C isotope ratio were able to model 
the significant decrease of the local ratio over the lifetime of the 
sun (about a factor two), which was thought at that time to be 
present. However, the models require some fine tuning and cali- 
bration, so that a smaller decrease is not in obvious conflict with 
the models. Nevertheless, some variation of the '^C/^^C isotope 
ratio is to be expected. New models of the chemical evolution 
of the Galaxy are required to take into account the new observa- 
tional results. 

A larger sample of stars will allow a better estimate of the 
average ISM value. In addition, these observations will allow to 
determine the scale of the variations of the isotope ratio and thus 
the efficiency of mixing on different scales. 
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